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Sustainable  development  is  closely  associated  with  the  use  of  renewable  energy  resources.  In  order  to 
achieve  a  viable  development,  from  an  environmental  point  of  view,  the  energy  efficiencies  of  processes 
can  be  increased  using  renewable  energy  resources.  There  is  also  a  correlation  between  exergy  and 
sustainable  development,  since  exergy  is  consumed  or  destroyed  due  to  irreversibilities.  The  solar 
chimney  has  been  highlighted  in  studies  of  using  solar  energy  to  generate  electric  power.  In  this  paper, 
the  energy  and  exergy  analyses  of  the  airflow  inside  a  solar  chimney  are  presented.  Using  experimental 
data  obtained  in  a  prototype,  the  first  and  second  laws  of  thermodynamics  were  used  to  estimate  the 
amounts  of  energy  and  exergy  lost  to  the  surroundings  and  the  exergetic  efficiency.  The  dead  state  was 
defined  using  two  different  reference  temperatures.  The  results  show  that  the  exergy  losses  were  lower 
and  the  efficiency  was  higher  for  the  lowest  ambient  temperature  used  as  the  dead  state  temperature, 
when  compared  to  the  instantaneous  ambient  temperature. 
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1.  Introduction 

Solar  energy  is  a  source  of  clean  and  renewable  energy,  which 
produces  neither  greenhouse  gases  nor  hazardous  waste  through 
use.  Solar  energy  technologies  include  heating,  photovoltaic,  thermal 
electricity  and  solar  architecture.  In  usual  thermal  applications,  the 
solar  energy  is  applied  to  create  steam  to  drive  an  electrical 
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generator.  A  relatively  new  device  that  uses  the  solar  energy  to 
generate  electricity  is  a  solar  chimney.  The  solar  chimney  combines 
the  technologies  of  solar  collectors,  chimneys  and  turbines.  The  solar 
radiation  reaches  and  passes  through  the  cover,  heating  the  soil 
beneath  the  cover,  which  in  turn  heats  the  air.  The  buoyancy  forces 
caused  by  temperature  gradients  in  the  air  promote  upward  airflow 
in  the  tower  and  through  the  centre  of  the  lid,  and  this  airflow  can  be 
used  to  drive  turbines  to  produce  electricity  [1-3].  The  concept  was 
developed  and  a  prototype  was  operated  during  the  1980s  by  Schaich 
[2  .  The  first  practical  application  of  a  solar  chimney  was  the 
development  of  a  prototype  that  generated  50  kW  of  power  built 
in  Manzanares,  Spain  [4],  in  1981/1982.  The  plant  operated  in  the 
period  from  1982  to  1989,  and  was  connected  to  the  local  power 
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supply  between  1986  and  1989  [3  .  The  operation  and  results  of 
operational  tests  of  the  prototype  are  described  by  Haaf  et  al.  [5]  and 
Haaf  [6],  No  complete  solar  power  plant  using  a  solar  chimney  has 
been  built  so  far,  but  a  project  of  200  MW  is  being  developed  in 
Australia  [7]. 

Economic  evaluations  based  on  measurements  taken  and  the 
operational  experience  gained  during  the  Manzanares  project 
showed  that  100  MW  plants  are  economically  competitive  to 
conventional  power  plants  [8].  To  evaluate  the  viability  of  new 
projects,  cost  models  for  large-scale  solar  chimney  power  plants 
were  developed.  Schlaich  [3]  presented  energy  production  costs 
and  compared  these  values  to  coal  and  combined  cycle  power 
plants  based  on  equal  and  common  methods.  Schlaich  et  al.  [9] 
presented  cost  values  for  the  solar  chimney  components  and 
levelised  electricity  costs  for  various  plants  for  fixed  economic 
parameters.  Bernardes  [10]  determined  the  costs  for  plants  of 
several  sizes  and  investigated  the  sensitivity  of  the  levelised 
electricity  costs  to  economic  parameters.  Fluri  et  al.  11]  presented 
a  more  detailed  cost  model,  including  a  model  for  the  power 
conversion  unit  and  compared  the  results  to  the  literature  models. 

Many  works  have  been  developed  on  the  study  of  solar 
chimneys.  Zhou  et  al.  [12]  presented  a  comprehensive  state  of 
the  art  review  of  research  and  development  of  solar  chimneys  for 
power  generation.  The  authors  reviewed  experimental  and  theo¬ 
retical  studies,  concerning  both  the  physical  and  economic  aspects 
of  developments,  and  concluded  that  solar  chimney  power  tech¬ 
nology  will  possibly  play  an  important  role  in  world  development 
in  the  future. 

Most  theoretical  works  about  solar  chimneys  found  in  litera¬ 
ture  deal  only  with  the  analysis  based  on  the  first  law  of 
thermodynamics.  In  recent  years,  growing  attention  has  been 
given  to  analysis  of  drying  processes  using  the  second  law  of 
thermodynamics  [13-19].  However,  the  only  research  developed 
with  a  second  law  of  thermodynamics  analysis  for  a  SCPP  was 
performed  by  Petela  [20].  The  author  presented  a  thermodynamic 
analysis  of  energy  distribution  in  the  whole  system  using  an 
exergetic  approach.  The  objective  was  to  outline  a  simplified 
interpretative  mathematical  model  of  the  solar  chimney.  However, 
the  results  presented  were  based  on  a  series  of  assumptions,  not 
on  experimental  data.  The  main  objective  of  the  present  paper  is 
to  analyse  the  thermodynamic  behaviour  of  a  solar  chimney 
prototype  without  any  products  inside.  Applying  the  first  and 
second  laws  of  thermodynamics,  the  energy  transferred  to  the 
airflow  and  the  exergy  lost  to  the  surroundings,  as  well  as  the 
exergetic  efficiency,  and  the  energy  and  exergy  rates  inside  the 
chimney  are  estimated.  The  results  obtained  and  presented  herein 
are  from  experimental  data  from  a  4-day  test  in  which  the  solar 
radiation,  ambient,  ground  and  outlet  airflow  temperatures,  inlet 
and  outlet  humidities  and  the  outlet  mass  flow  rate  were  mea¬ 
sured.  It  is  important  to  mention  that  the  prototype  has  small 
dimensions  in  comparison  to  the  dimensions  of  solar  chimney 
power  plants,  since  the  main  purpose  of  the  prototype  was  the 
drying  of  agricultural  products. 


2.  Literature  review 

There  is  now  a  wide  and  growing  body  of  literature  on  the 
study  of  solar  chimneys.  The  outstanding  works  are  discussed  in 
this  section. 

2.1.  Mathematical  models 

Several  studies  about  solar  chimneys  have  been  performed  to 
predict  the  behaviour  of  the  airflow  inside  a  solar  chimney,  using 


both  computational  fluid  dynamics  (CFD)  simulations  and  global 
analysis. 

Bernardes  et  al.  [21]  presented  a  theoretical  analysis  of  the 
laminar  natural  convection  inside  the  device,  for  steady  state 
conditions.  After  this  study,  Maia  et  al.  [22  improved  this  model 
to  take  into  account  unsteady  turbulent  conditions,  validating  the 
results  through  comparison  with  experimental  data.  Chergui  et  al. 
[23]  developed  a  finite  volume  method  to  discretise  the  Navier- 
Stokes  and  energy  equations  for  a  laminar  natural  convection 
airflow.  A  mesh  size  of  1/22  (24x24  cells)  was  used  and  the 
results  were  validated  with  respect  to  the  Vahl  Davis  benchmark 
solution  [24].  Pastohr  et  al.  [25]  modelled  the  problem  using  the 
software  Fluent.  The  Reynolds-averaged  Navier-Stokes  equations 
were  numerically  solved  for  the  airflow,  with  a  simple  turbine 
model  and  the  heat  conduction  equation  to  predict  the  ground 
surface  temperature.  Ming  et  al.  [26]  solved  the  Navier-Stokes  and 
energy  equations  for  the  airflow,  as  well  as  transport  equations  for 
the  turbulent  variables.  Different  mathematical  models  for  the 
collector,  the  chimney,  and  the  energy  storage  layer  under  the 
solar  chimney  were  developed  to  be  used  as  boundary  conditions 
for  the  aforementioned  equations.  The  Manzanares  prototype  was 
selected  as  a  physical  model  for  the  numerical  simulation.  Tingz- 
hen  et  al.  [27]  carried  out  numerical  simulations  on  solar  chimney 
systems  using  the  CFD  techniques.  The  influence  of  solar  radiation 
and  pressure  drop  across  the  turbine  on  the  flow  and  heat  transfer, 
output  power,  and  energy  loss  of  the  solar  chimney  power  plant 
system  (SCPP)  were  analysed.  The  numerical  simulation  results 
revealed  that,  for  a  solar  chimney  with  a  chimney  400  m  in  height 
and  with  a  collector  1500  m  in  radius  and  a  5-blade  turbine,  the 
maximum  power  output  and  turbine  efficiency  is  about  10  MW 
and  50%,  respectively. 

Koonsrisuk  and  Chitsomboon  28]  showed  that  dynamic  simi¬ 
larity  for  scaled  models  and  a  full-scale  solar  chimney  prototype  is 
achievable,  while  maintaining  the  same  insolation.  Nevertheless, 
only  a  partial  geometrical  similarity  is  allowed.  Experiments  were 
not  performed  to  confirm  the  validity  of  the  study.  Instead,  data 
similarities  between  the  model  and  a  prototype  were  compared 
using  numerical  results  provided  by  CFD  analysis.  Later,  the  same 
authors  [29]  used  dimensional  analysis  to  combine  eight  primitive 
variables  into  one  dimensionless  parameter  to  establish  dynamic 
similarity  between  a  prototype  and  its  scaled  models.  This  dimen¬ 
sionless  parameter  was  interpreted  as  the  total  kinetic  energy 
scaled  by  the  buoyant  energy  of  the  rising  hot  air.  To  validate  the 
proposed  dimensionless  parameter,  three  classes  of  scaled  models 
were  set  up  and  the  governing  equations  were  solved  using  CFD 
analysis.  Xu  et  al.  [1]  developed  a  two  dimensional  model  for  the 
solar  chimney  and  an  energy  storage  layer  under  the  device, 
evaluating  the  influences  of  solar  radiation,  turbine  pressure  drop 
and  efficiency  on  the  flow  and  heat  transfer,  output  power,  and 
energy  loss  of  the  system.  The  Fluent  software  was  used  to  solve 
the  governing  equations  and  the  model  was  validated  with 
experimental  results  from  the  Manzanares  prototype.  Sangi  et  al. 
[30]  described  a  detailed  numerical  analysis  of  a  SCPP,  using  Fluent 
software  to  simulate  a  two-dimensional  axisymmetric  model,  with 
the  standard  k-epsilon  turbulence  model.  The  model  was  validated 
with  data  from  the  Manzanares  prototype.  The  authors  also 
developed  an  expression  for  the  static  pressure  inside  the 
collector. 

The  airflow  generated  in  a  solar  chimney  is  commonly  used  to 
drive  a  wind  turbine  and  to  generate  electricity.  Therefore,  the 
electrical  output  is  the  most  important  characteristic  in  a  SCPP. 
The  prediction  of  its  performance  and  the  influence  of  the  physical 
dimensions  on  a  solar  chimney's  performance  have  been  studied 
by  several  researchers.  Asnaghi  and  Ladjevardi  [31]  stated  that  the 
most  important  physical  elements  in  a  solar  chimney  are  the 
tower  dimensions  as  they  cause  the  most  significant  influence  on 
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the  flow  patterns.  Another  important  factor  is  the  collector's  roof 
shape.  Also,  research  has  shown  that  greater  power  production  is 
feasible  by  optimising  the  collector's  roof  shape  and  height. 
Bernardes  et  al.  [32]  developed  a  model  to  estimate  the  power 
output  of  solar  chimneys  and  examine  the  effects  of  various 
ambient  conditions  and  structural  dimensions  on  the  power  out¬ 
put.  The  model  was  validated  using  results  from  the  Manzanares 
pilot  plant.  Koonsrisuk  and  Chitsomboon  [33]  compared  the 
predictions  of  performances  of  SCPP  by  using  five  simple  theore¬ 
tical  models  proposed  in  the  literature.  The  authors  performed 
CFD  simulations  and  compared  the  results  with  theoretical  pre¬ 
dictions.  Analytical  approaches  are  also  found  in  literature  for 
estimation  of  the  airflow  parameters  of  SCPP.  Pretorius  and  Kroger 
[34]  evaluated  the  influence  of  several  convective  heat  transfer 
coefficients,  turbine  inlet  loss  coefficients,  collector  roof  glasses 
and  various  types  of  soil  on  the  performance  of  a  SCPP.  Fluri  and 
Von  Backstrom  [35]  determined  the  SCPP  power  production 
spectrum  over  a  year  and  determined  the  required  characteristics 
of  the  turbine(s)  to  match  the  plant  and  compared  three  power 
conversion  unit  configurations.  Nizetic  and  Klarin  [36]  presented  a 
mathematical  model  based  on  a  simplified  thermodynamic  ana¬ 
lysis  of  the  overall  solar  chimney  cycle  to  estimate  the  optimal 
factor  of  turbine  pressure  drop.  The  authors  theorised  that  even 
for  constant  solar  irradiation,  the  pressure  potential  is  not  fixed, 
but  a  function  of  the  air  temperature  increase  in  the  collector. 
Li  et  al.  [37]  described  a  comprehensive  theoretical  model  for  the 
performance  evaluation  of  a  solar  chimney  power  plant,  verified 
by  the  experimental  data  from  the  Spanish  prototype.  The  authors 
stated  that  there  is  a  maximum  power  output  under  a  given  solar 
radiation  condition  and  that  there  is  a  limitation  on  the  maximum 
collector  radius  for  the  maximum  attainable  power  of  the  SCPP; 
whereas  no  such  limitation  exists  for  chimney  height  in  terms  of 
contemporary  construction  technology.  Hamdan  [38;  provided  a 
mathematical  thermal  model  to  evaluate  the  performance  of  a 
SCPP,  based  on  its  dimensions  and  operating  conditions.  The 
model  was  developed  using  a  modified  Bernoulli  equation  with 
buoyancy  effect  and  an  ideal  gas  equation.  Koonsrisuk  and 
Chitsomboon  [39]  developed  a  theoretical  model  to  evaluate  the 
performance  of  SCPP.  The  relationships  between  the  pressure  ratio 
and  the  mass  flow  rate  and  between  the  temperature  rise  across 
the  collector  and  the  power  output  were  also  presented. 

2.2.  Worldwide  studies 

Many  attempts  have  been  made  to  evaluate  solar  chimneys  in 
some  parts  of  the  world.  To  the  best  knowledge  of  the  authors, 
there  was  only  one  prototype  built  in  South  America.  Maia  et  al. 
[40]  described  an  experimental  study  of  a  prototype  built  in  Brazil 
with  a  solar  collector  with  a  25  m  diameter  and  a  13  m  high  tower. 
Since  the  estimated  power  was  very  low,  a  turbine  was  not 
installed.  Instead,  the  hot  airflow  generated  was  used  to  dry 
agricultural  products.  Ferreira  et  al.  [41  investigated  the  technical 
feasibility  of  the  device  as  a  dryer.  Drying  tests  of  coffee  grains, 
bananas  and  tomatoes  were  performed  inside  the  dryer. 

Kasaeian  et  al.  [42]  conducted  an  experimental  investigation  of 
climatic  effects  on  the  efficiency  of  a  SCPP.  The  solar  chimney  in 
pilot  scale  was  built  in  Iran,  with  a  chimney  height  of  12  m  and  a 
collector  with  a  10  m  diameter.  Sangi  [43]  estimated  the  potential 
for  electrical  energy  production  by  a  solar  chimney  in  selected 
regions  of  Iran.  A  simplified  mathematical  model  based  on  the 
energy  balance  was  developed  to  estimate  the  power  output  of 
solar  chimneys,  as  well  as  to  examine  the  effect  of  various  ambient 
conditions  and  structural  dimensions  on  power  generation.  Asna- 
ghi  and  Ladjevardi  [31]  performed  numerical  simulations  to 
investigate  SCPP  in  12  selected  locations  in  Iran,  comparing  power 
generation  for  each  region.  An  analysis  of  overall  performance  of  a 


device  with  similar  fundamental  dimensions  to  the  Manzanares 
prototype  was  provided,  considering  the  radiation  variations 
through  the  year  in  the  selected  regions.  Gholamalizadeh  and 
Mansouri  44]  presented  a  comprehensive  analysis  to  predict  the 
performance  of  a  SCPP  in  Kerman,  Iran.  Analytical  and  numerical 
models  were  developed,  and  the  results  were  compared  to 
experimental  data  from  the  Manzanares  prototype. 

Dai  et  al.  [2]  evaluated  a  SCPP  to  be  constructed  in  north¬ 
western  China  to  provide  electric  power  for  remote  villages. 
The  performance  of  the  device  was  compared  for  three  possible 
locations.  Results  showed  that  a  prototype  200  m  high  and  with  a 
500  m  collector  diameter  was  expected  to  produce  110-190  kW 
electric  power.  Zhou  et  al.  [45]  developed  a  small-scale  prototype 
of  solar  chimney  thermal  power  plant  in  China,  8  m  high  and  with 
a  radius  that  may  vary  from  1  m  to  5  m.  A  simulation  study  was 
carried  out  to  investigate  the  performance  of  the  system  based  on 
a  mathematical  model.  With  the  reduced  dimensions,  the  device 
was  able  to  generate  up  to  approximately  5.5  W.  Zhou  et  al.  [46] 
built  a  prototype  in  China,  consisting  of  an  air  collector  10  m  in 
diameter  and  an  8  m  tall  chimney.  The  temperature  distribution 
was  presented  for  a  typical  cool  day  and  for  a  typical  warm  day. 
The  results  showed  a  temperature  difference  between  the  collec¬ 
tor  outlet  and  the  ambient  up  to  24.1  °C.  Zhou  et  al.  [47]  developed 
a  mathematical  model  for  the  evaluation  of  the  performance  of  a 
solar  chimney  to  generate  electric  power  to  run  the  railway  lines 
in  the  Qinghai-Tibet  Plateau  and  to  supply  the  developed  regions 
by  high-voltage  lines.  The  results  indicated  that  when  10-20%  of 
land  in  the  plateau  is  used  for  solar  chimney  power  plants,  the 
yearly  power  could  account  for  10.7-21.3%  of  China's  energy 
consumption  in  2008. 

Onyango  and  Ochieng  [48]  evaluated  the  potential  of  solar 
chimneys  for  application  in  developing  countries,  highlighting  its 
appropriateness  to  rural  villages.  Nizeticet  al  [49]  evaluated  the 
viability  of  implementing  SCPP  in  the  Mediterranean  region 
through  a  cost  analysis.  A  device  550  m  high  with  a  collector  with 
a  1250  m  diameter  would  generate  2.8-6.2  MW  of  power.  Larbi 
et  al.  [50]  presented  a  performance  analysis  of  a  prototype 
expected  to  produce  between  140  kW  and  200  kW  of  electricity 
in  Algeria's  south-western  region.  Hamdan  [7]  developed  a  sim¬ 
plified  thermodynamic  analytical  model  to  analyse  the  feasibility 
of  solar  chimney  power  plants  for  the  Arabian  Gulf  region. 
Ketlogetswe  et  al.  [51]  described  an  experimental  study  of  a  solar 
chimney  with  a  collection  area  of  approximately  160  m2  built  in 
Botswana  (Southern  Africa).  The  tower  had  an  inner  diameter  of 
2  m  and  a  height  of  22  m.  The  main  parameters  evaluated  were  the 
air  velocity,  temperature  and  solar  radiation.  The  results  were 
obtained  under  different  operational  conditions:  with  a  turbine 
installed,  without  the  turbine  and  without  the  turbine  and  the 
diffuser.  The  range  of  temperature  difference  varied  from  about  2.1 
to  7.5  °C.  Air  velocity  was  in  the  range  of  1-2.5  m/s  with  the 
diffuser  installed  and  2-4  m/s  with  the  diffuser  removed.  Other 
small-scale  prototypes  were  developed.  Kirsst  [52  built  a  10  W 
prototype  in  America.  Kulunk  [53]  built  a  2  m  high  prototype  in 
Turkey,  able  to  generate  0.14  W. 


3.  Experimental  setup 

A  prototype  of  a  solar  chimney  was  built  in  Belo  Horizonte 
(Brazil).  The  tower  was  12.3  m  high,  and  the  collector  diameter 
was  25  m.  The  height  of  the  collector  varied  from  0.05  m  in  the 
entrance  to  0.5  m  in  the  symmetry  axis  (Fig.  1). 

The  temperature,  velocity  and  humidity  ratio  of  the  air,  and  the 
ambient  temperature  and  solar  radiation  were  measured  over  a  4- 
day  period.  Despite  several  papers  (e.g.,  42,46])  presenting  results 
for  a  1-day  period  or  instantaneous  data,  it  was  chosen  to  show 
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the  results  for  a  4-day  period  in  order  to  highlight  the  influence  of 
the  climatic  parameters  on  the  airflow.  The  period  chosen  shows 
opposite  behaviours  concerning  the  incident  solar  radiation  and 
clearness  indexes.  The  ambient,  ground  and  outlet  airflow  tem¬ 
peratures  were  measured  with  I<  thermocouples  (uncertainty  of 
1  °C).  Eppley  Black  and  White  Pyranometers  Model  8-48  were 
used  to  determine  solar  radiation  (global  uncertainty  of  5%). 
Propeller  Homis  anemometers  were  used  to  measure  the  velocity 
in  the  tower.  The  global  uncertainty  was  6%.  Capacitive  psychrom- 
eters  were  used  to  measure  the  ambient  relative  humidity  and  the 
relative  humidity  of  the  airflow  (uncertainty  of  6%).  The  analogue 
voltage  or  current  signals  from  the  humidity,  temperature,  velocity 
and  solar  radiation  measurement  sensors  were  converted  into 
digital  signals  through  ADAMS  4018  Modules  in  a  data  acquisition 
system. 

The  ambient  temperature  and  relative  humidity  were  mea¬ 
sured  outside  the  device.  Thermocouples  and  psychrometers  were 
also  positioned  in  the  inlet  and  outlet  of  the  device.  Anemometers 
were  used  to  measure  the  velocity  along  a  cross  section  of  the 
tower.  A  numerical  integration  was  performed  to  determine  the 
mass  flow  rate  of  the  airflow.  The  global  and  diffuse  components 
of  solar  radiation  were  measured  in  a  horizontal  plane  outside  the 
prototype. 

The  experiments  were  conducted  in  late  autumn  in  Brazil,  with 
maximum  and  minimum  air  temperatures  around  27.3  °C  and  8.4  °C, 
respectively.  The  air's  relative  humidity  varied  between  25.6%  and 
90.6%  and  the  maximum  solar  radiation  was  1065  W/m2. 


4.  Modelling 

The  balance  equations  set  were  developed  according  to  models 
from  the  literature  15,54].  The  balance  equations  were  applied  to 
find  the  work  and  heat  interactions,  the  exergy  rates,  the  exergy 
losses,  and  the  exergy  efficiency. 

The  general  equation  of  mass  conservation  can  be  expressed  in 
the  rate  form  as 

mai  =  mao  0) 

where  mai  and  mao  represent  the  inlet  and  outlet  mass  flows  of  the 
airflow,  respectively. 

For  the  drying  air  and  moisture,  the  mass  conservation  is 
expressed  as 

mWi  +  rhmp  =  mwo  (2) 

mwi  and  rhwo  represent  the  inlet  and  outlet  mass  flows  of  humidity, 
respectively.  mmp  is  the  mass  flow  of  the  moisture  from  the 
ground.  The  mass  conservation  for  the  moisture  can  also  be 


rewritten  in  terms  of  the  inflow  and  outflow  specific  humidities, 
respectively,  and  mao ,  as 


bblai(J)ai  "T  Wimp  —  WlaoCOao 


(3) 


The  general  energy  balance  can  be  expressed  as  the  total 
energy  inputs  equal  to  total  energy  outputs: 


l£fn  = 


out 


(4) 


or 


v: 


Q, =  wiao  (  ha0  -f  dp-  ] -mai  hai  + 


(5) 


The  previous  equation  is  used  to  determine  the  net  heat  rate  Q. 
The  energy  utilisation  rate  is  neglected,  since  there  is  no  turbine  in 
the  system.  Vat  and  Vao  represent  the  air  velocity  in  the  inlet  and 
outlet  of  the  system,  respectively.  The  specific  enthalpies  of  the  air 
in  the  inlet  and  outlet,  hai  and  hao,  were  determined  using  the 
correlations  available  in  the  EES  (Engineering  Equation  Solver) 
software: 

The  general  exergy  balance  can  be  written  as  follows: 


(6) 


In  the  previous  equation,  Exin,  Exout and  Exiost  represent  the 
exergy  inflow,  outflow,  and  the  exergy  loss  respectively.  This 
equation  can  be  also  expressed  as 

EXiieaf—Exwork  +  Exmassjn— EXmass, out  —  EX[0S  f  (7) 

With  the  exergy  flow  rate  due  to  the  heat  transfer  given  as 


o 


EXheat  =  (1-^  IO¬ 


CS) 


Q  is  the  heat  transfer  rate  through  the  boundary  at  temperature  Tk 
at  location  k,  considered  as  the  ground  temperature.  T0  represents 
the  dead  state  temperature.  For  the  system  considered,  the  exergy 
rate  due  to  work  interactions  Exwork  can  be  neglected. 

The  exergy  inflow  Exmassin  is  only  due  to  the  airflow  entering 
the  system,  and  the  exergy  outflow  Exmass,out  is  due  to  the  airflow 
leaving  the  system  and  to  the  water  removed  from  the  ground, 
Exw. 


Exmass[n  —  WlaiWai 

Exmass, out  —  EflaoW  ao  +  Exw 

where 


(9) 


(10) 


(11) 


EXW  —  WlmpWwo 

Therefore,  the  exergy  outflow  is 

EXmass,  out  —  WlaoW ao  +  TTlmpy/w0  (12) 

The  specific  flow  exergy  *F  is  calculated  in  a  generic  form  as 

V  =  (h-ho)-T0(s-s0)  (13) 

where  h  is  the  specific  enthalpy  and  s  is  the  specific  entropy. 

The  inlet  and  outlet  specific  flow  exergies  are  calculated  by  55] 


Wai  —  (Ep^ai  +  toaiCpy)T  q 


Tai  -l-ln^i)  +(1  +  1.6078<wa,)^T0ln^ 


0 


0 


0 


+Ra  To 


(1  +  1.6078„„,)ln| 


COai 

COq 


(14) 


Wao  ~  ( Ep,ao  +  <WaoEp^)To 


1-ln 


+(1  +1.6078fflao)KflT0ln^ 

I  0 
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+^a  T  o 


(1  +  1.6078©ao)ln 


n  +  1.6078©0\ 
VI  +  1.6078©aoJ 


4-  1.6078®aoln  — 

COo  _ 

(15) 


where  Ra  is  the  ideal  air  constant,  PQ  is  the  dead  state  pressure,  Pai 
and  Pao  are  the  air  pressure  in  the  inlet  and  outlet,  respectively, 
and  Tai  and  Tao  are  the  air  temperatures  in  the  inlet  and  outlet, 
respectively.  Cpy  is  the  specific  heat  of  the  water  vapour,  and  &>0  is 
the  specific  humidity  of  the  flow  in  dead  state. 

The  outlet  specific  flow  exergy  of  the  water,  considered  as  an 
incompressible  substance,  is  given  by 

Vwo  =  c(rao-T0-To\n1£\  (16) 


It  is  assumed  that  the  air  and  water  leave  the  system  at  the 
same  temperature,  Tao.  In  the  previous  equation,  C  is  the  specific 
heat  of  water. 

Therefore,  the  exergy  loss  is  determined  by 

EX/ost  =  Q  ^aiWai  th aoWao  ^mpWwo  07) 

The  irreversibility  rate  /  is  given  by 

I  =  Exl0St  (18) 


The  exergy  efficiency  e  is  defined  as  the  ratio  of  total  exergy 
output  to  total  exergy  input: 


Exout 
e  =  — - 

EXin 


=  1 


Exiost 

EXin 


(19) 


Van  Gool  [56]  proposed  that  maximum  improvement  in  the 
exergy  efficiency  for  a  process  or  system  is  achieved  when  the 
exergy  loss  or  irreversibility  is  minimised.  Therefore,  it  is  useful  to 
use  the  concept  of  exergetic  improvement  potential,  defined  as 


IP  —  (1 —  £)(EXin—  Exout) 


(20) 


Exergy  is  always  evaluated  with  respect  to  a  reference  envir¬ 
onment,  the  dead  state.  When  a  system  is  in  thermodynamic 
equilibrium  with  the  environment,  the  state  of  the  system  is  called 
the  dead  state  due  to  the  fact  that  the  exergy  is  zero.  According  to 
[54],  dead  state  is  an  arbitrary  reference  state.  Usually,  results  are 
presented  in  literature  for  steady  state  conditions.  In  this  paper, 
unsteady  state  conditions  were  evaluated.  The  reference  tempera¬ 
ture  (selected  as  the  dead  state  temperature)  was  the  inlet 
temperature.  Since  the  inlet  temperature  varied  during  the  tests, 
two  dead  states  were  defined  and  the  results  were  compared.  Both 
dead  states  were  selected  as  a  system  with  zero  velocity  and 
elevation  relative  to  coordinates  in  the  environment,  with  the  local 
ambient  atmospheric  pressure  (91.5  kPa).  The  main  difference 
between  the  definitions  was  the  temperature:  in  the  first,  the 
selected  temperature  was  the  minimum  ambient  temperature 
(approximately  8.5  °C)  during  the  experimental  tests  performed; 
in  the  second  dead  state,  the  temperature  was  the  instantaneous 
ambient  temperature.  It  can  be  noticed  that  in  all  the  above 
balance  equations  (Eqs.  (2)— (6))  the  cumulative  terms  are 
neglected,  considering  that  even  if  the  thermodynamics  properties 
change  with  time,  its  temporal  derivatives  are  null. 


5.  Results  and  discussion 

Fig.  2  presents  the  total  and  diffuse  components  of  solar 
radiation  measured  during  the  four  days  of  tests.  The  average 
clearness  index  (the  ratio  of  the  daily  radiation  to  the  extra¬ 
terrestrial  radiation  of  the  day  [57  )  was  determined  for  each  day, 
obtaining  values  of  0.68,  0.63,  0.72  and  0.53,  respectively,  for  the 
first,  second,  third  and  fourth  day  of  tests.  It  can  be  seen  that 
higher  clearness  indexes  correspond  to  clearer  days. 

Fig.  3  shows  the  inlet  and  outlet  temperatures  of  the  airflow 
and  the  ground  temperature.  As  expected,  the  temperature  and 
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Fig.  2.  Total  and  diffuse  components  of  solar  radiation. 
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Fig.  3.  Temperature  distribution. 
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Fig.  4.  Specific  humidity  distribution. 


total  solar  radiation  profiles  exhibit  a  similar  behaviour.  Higher 
temperatures  are  observed  for  higher  clearness  indexes.  It  can  also 
be  seen  that  the  solar  chimney  is  able  to  operate  at  night.  During 
the  day,  when  there  is  incidence  of  solar  radiation,  ground  surface 
temperature  is  higher  than  the  temperature  of  the  deeper  layers. 
Therefore,  heat  is  conducted  from  the  surface  to  the  deeper  layers  of 
the  ground  and  stored  there  as  heat.  At  night,  the  temperature  of  the 
ground  surface  decreases,  becoming  lower  than  the  deeper  layers' 
temperatures.  Thus,  heat  flux  goes  in  the  opposite  direction  [40]. 

Fig.  4  shows  the  inlet  and  outlet  airflow  specific  humidities. 
It  can  be  seen  that  the  outlet  humidity  is  greater  than  the  inlet 
humidity.  Since  the  prototype  ground  was  not  sealed,  the  differ¬ 
ence  can  be  explained  by  the  removal  of  water  from  the  ground  by 
the  airflow.  A  difference  between  the  inlet  and  outlet  airflow 
specific  humidities  is  also  observed.  During  the  day,  the  ambient 
temperatures  were  higher,  causing  a  decrease  in  the  specific 
humidity  in  the  inlet.  Since  the  air  entered  the  device  with  a 
lower  specific  humidity,  it  was  easier  to  remove  water  from  the 
ground,  increasing  the  difference  in  the  specific  humidities  in  the 
inlet  and  outlet.  At  night,  the  opposite  behaviour  was  found. 

The  airflow  is  produced  by  buoyancy  forces  resulting  from  a 
greenhouse  effect  inside  the  collector.  Therefore,  the  velocity 


depends  on  the  difference  between  the  ground  and  the  airflow 
temperature,  which  depends  on  the  incident  solar  radiation.  Fig.  5 
presents  the  mass  flow  rate  during  the  whole  test.  The  higher 
values  were  obtained  when  the  solar  radiation  and  the  tempera¬ 
ture  were  higher,  but  there  is  also  an  airflow  at  night,  due  to  the 
temperature  difference  observed  in  Fig.  3.  The  average  velocity  in 
the  tower  varied  from  1.1  m/s  to  2.7  m/s,  corresponding  to 
Reynolds  numbers  (in  the  tower)  from  6.4  x  104  to  1.5  x  105, 
showing  characteristics  of  turbulent  airflows. 

Fig.  6  shows  the  heat  transfer  rate  from  the  ground  to  the 
airflow.  As  expected,  the  heat  transfer  rate  exhibits  a  similar 
behaviour  to  the  other  parameters,  with  higher  values  associated 
to  higher  solar  radiation  values.  There  is  also  a  small  portion  of 
heat  being  transferred  at  night,  when  there  is  no  incidence  of  solar 
radiation. 

As  previously  discussed,  the  dead  state  is  an  arbitrary  state  of 
reference.  In  this  work,  two  dead  states  were  selected.  Fig.  7 
presents  the  exergy  rates  for  a  fixed  value  of  the  temperature  of 
the  dead  state.  The  temperature  was  selected  as  the  minimum 
value  of  the  ambient  temperature  obtained  during  the  tests.  Fig.  8 
presents  the  exergy  rates  for  the  dead  state  temperature  defined 
as  the  ambient  temperature,  which  varied  from  8.4  °C  to  27.3  °C. 


Exergy  rates  (W)  Heat  transfer  rate  (kW)  Mass  flow  rate  (kg/s) 
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Fig.  5.  Mass  flow  rate  distribution. 
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Fig.  6.  Heat  transfer  rate  distribution. 
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Fig.  7.  Exergy  as  a  function  of  time  for  the  dead  state  selected  with  a  fixed  temperature. 
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Fig.  8.  Exergy  as  a  function  of  time  for  the  dead  state  selected  with  ambient  temperature. 
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Fig.  9.  Exergy  losses  as  a  function  of  time  for  the  dead  state  selected  with  a  fixed  temperature. 


The  general  behaviour  of  the  exergy  rates  is  similar  to  the 
behaviour  of  the  temperature.  The  exergy  rates  are  higher  for  a 
constant  temperature  of  the  dead  state,  since  a  higher  amount  of 
work  could  be  produced  as  the  system  comes  to  equilibrium  with 
the  reference  environment  at  a  lower  temperature.  In  Figs.  7  and  8, 
the  higher  exergy  rates  were  found  for  the  exergy  flow  rates  due  to 
heat  transfer,  and  the  lower  exergy  rates  were  found  for  the  exergy 
rates  of  the  water  removed  from  the  ground.  Therefore,  the  higher 
exergy  rates  are  associated  to  higher  solar  radiation  levels.  When 
the  ambient  temperature  was  selected  as  the  dead  state  tempera¬ 
ture,  the  exergy  rate  of  the  inlet  mass  flow  is  always  zero,  since  the 
reference  state  is  the  inlet  air. 

Figs.  9  and  10  present  the  exergy  losses  for  a  fixed  value  of  the 
temperature  of  the  dead  state  and  for  the  dead  state  temperature 
defined  as  the  ambient  temperature,  respectively.  The  losses  are 
very  significant,  and  the  higher  losses  correspond  to  higher  solar 
radiation  levels,  due  to  the  higher  heat  transfer  rates  achieved. 
Since  the  hot  airflow  generated  was  only  used  to  remove  humidity 
from  the  ground,  not  being  used  to  drive  any  engine  or  to  dry 
agricultural  products,  the  heat  was  not  used,  increasing  the 
exergetic  losses.  It  is  also  seen  that  the  losses  are  more 


pronounced  for  a  fixed  dead  state  temperature  due  to  the  lower 
comparison  level  as  it  is  more  difficult  to  enter  into  equilibrium 
with  a  lower  reference  environment  and/or  because  the  exergy 
rates  (mainly  related  to  the  exergy  rates  due  to  heat  transfer)  are 
higher  in  this  case. 

Figs.  11  and  12  present  the  exergetic  efficiency  for  a  fixed  value 
of  the  temperature  of  the  dead  state  and  for  the  dead  state 
temperature  defined  as  the  ambient  temperature,  respectively.  It 
can  be  seen  that  the  efficiencies  are  higher  for  a  fixed  value  of  the 
dead  state  temperature,  since  the  reference  temperature  is  lower 
than  the  ambient  temperature  for  the  majority  of  the  time.  When 
analysing  the  temporal  change  in  temperatures,  it  can  be  seen  that 
the  exergetic  efficiency  is  continuously  varying,  with  higher  values 
obtained  for  the  last  day  of  tests.  In  the  last  day,  lower  solar 
radiation  levels  were  found,  as  well  as  lower  heat  transfer  rates. 
Nevertheless,  the  exergetic  losses  were  lower,  increasing  the 
exergetic  efficiency. 

No  data  in  the  literature  was  found  for  a  similar  configuration 
of  the  prototype  used  in  these  experiments.  Hamdan  [7]  presented 
a  simplified  thermodynamics  analytical  model  for  steady  airflow 
inside  a  solar  chimney  power  plant.  In  order  to  validate  the 
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Fig.  10.  Exergy  losses  as  a  function  of  time  for  the  dead  state  selected  with  ambient  temperature. 
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Fig.  11.  Exergetic  efficiency  as  a  function  of  time  for  the  dead  state  selected  with  a  fixed  temperature. 


developed  model,  the  results  were  compared  against  published 
experimental  data,  with  good  agreement  for  the  power  generated 
as  a  function  of  geometric  and  operational  parameters.  The 
second-law  efficiency  of  the  system  was  determined,  based  on 
the  developed  model,  with  values  between  15%  and  38%,  depend¬ 
ing  on  the  chimney  height  (between  500  m  and  1000  m).  It  was 
observed  that  the  second-law  efficiency  increases  with  the  chim¬ 
ney  height.  The  efficiency  defined  is  related  to  the  expected 
turbine  power,  as  also  suggested  by  Saidur  et  al.  [58]  for  solar 
power  plants.  Nevertheless,  the  results  obtained  by  Hamdan  [7] 
cannot  be  compared  with  the  efficiencies  obtained  in  this  paper, 
since,  in  this  paper,  the  solar  chimney  does  not  have  a  turbine  and 
the  dimensions  of  the  prototypes  were  not  alike. 

Figs.  13  and  14  present  the  improvement  potential  for  a  fixed 
value  of  the  temperature  of  the  dead  state  and  for  the  dead  state 
temperature  defined  as  the  ambient  temperature,  respectively.  As 
expected,  the  improvement  potential  is  lower  for  a  fixed  tempera¬ 
ture  of  the  dead  state,  since  the  efficiencies  are  higher  for  this 
condition.  The  relatively  higher  values  of  the  improvement 


potential  can  be  explained  by  the  fact  that  the  airflow  is  only 
being  used  to  remove  water  from  the  ground.  Therefore,  there  is  a 
great  potential  of  exergy  to  be  harnessed,  reducing  the  exergy 
losses. 

As  a  core  outcome  of  the  results,  it  can  be  said  that  higher  solar 
radiation  leads  to  higher  temperatures  and  mass  flow  rates, 
resulting  in  higher  the  heat  transfer  rates.  Since  the  higher  exergy 
rates  are  due  to  the  heat  transfer  rates,  the  higher  exergy  losses 
and  the  lower  exergetic  efficiencies  can  be  related  to  the  higher 
solar  radiation  levels.  However,  it  must  be  clear  that  these  results 
are  valid  for  this  prototype,  where  the  heat  transfer  rate  was  only 
used  to  remove  water  from  the  ground. 

The  temperature  defined  for  the  dead  state  played  an  impor¬ 
tant  role  on  the  exergetic  analysis.  When  a  fixed  temperature  was 
selected,  instead  of  the  ambient  temperature  (which  varied  during 
the  tests),  higher  exergy  rates  and  exergetic  efficiencies  were 
found.  The  exergetic  losses  were  significant  due  to  the  fact  that 
the  hot  airflow  was  only  used  to  remove  water  from  the  ground; 
nevertheless,  the  main  goal  of  this  paper  was  to  analyse  the 
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Fig.  12.  Exergetic  efficiency  as  a  function  of  time  for  the  dead  state  selected  with  ambient  temperature. 
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Fig.  13.  Improvement  potential  as  a  function  of  time  for  the  dead  state  selected  with  a  fixed  temperature. 


behaviour  of  the  airflow  inside  the  solar  chimney  with  no  products 
inside.  It  is  important  to  note  that,  in  spite  of  this,  the  hot  airflow 
generated  in  a  solar  chimney  can  be  used  to  drive  a  turbine  in  solar 
chimney  power  plants  or  to  dry  agricultural  products  in  a  smaller 
prototype.  In  a  solar  chimney  power  plant,  it  is  shown  that  higher 
solar  radiation  levels  increase  the  power  output  and  the  energetic 
efficiency  of  the  system  [12,28].  In  a  drying  process,  the  most 
important  parameter  of  the  drying  efficiency  is  the  difference 
between  the  relative  humidities  of  the  airflow  and  the  products  to 
be  dried.  An  increase  of  the  source  of  energy  (in  this  case,  the  solar 
energy)  increases  the  airflow  temperature,  which  decreases  the 
relative  humidity  of  the  air,  improving  the  drying  characteristics  [59]. 


6.  Conclusions 

In  order  to  perform  an  exergetic  and  energetic  analysis  of  the 
airflow  in  a  solar  chimney,  experimental  data  from  a  12.3  m  high 


prototype  of  a  solar  chimney  was  used.  The  tests  were  carried  out 
over  a  4-day  period.  The  main  points  of  the  conclusion  are 

“  the  profiles  of  temperature  and  velocity  are  closely  related  to 
incident  solar  radiation  on  the  device; 

-  the  temperature  defined  for  the  dead  state  plays  an  important 
role  in  the  exergy  rates  and  exergetic  efficiencies  obtained.  For 
a  fixed  temperature  of  the  dead  state  (defined  as  the  minimum 
ambient  temperature  of  the  whole  test),  it  was  found  that  the 
exergy  rates  and  exergetic  efficiencies  were  higher  than  those 
obtained  with  the  instantaneous  ambient  temperature.  This 
behaviour  can  be  explained  due  to  the  lower  reference 
environment; 

-  the  higher  exergy  rates  are  found  for  the  exergy  flow  rates  due 
to  heat  transfer; 

-  lower  exergy  losses  were  found  for  the  last  day.  Besides,  the 
higher  exergy  losses  were  found  close  to  12:00.  Therefore,  it 
can  be  concluded  that  higher  incident  solar  radiation  levels 
lead  to  higher  exergy  losses.  Higher  solar  radiation  levels 
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Fig.  14.  Improvement  potential  as  a  function  of  time  for  the  dead  state  selected  with  ambient  temperature. 


generate  higher  heat  transfer  rates  to  the  airflow.  Since  this 
heat  was  not  used,  the  exergy  losses  were  higher,  and  the 
exergetic  efficiencies  were  lower.  In  a  solar  chimney  power 
plant  or  in  a  small  prototype  of  a  solar  chimney  used  to  dry 
agricultural  products,  higher  solar  radiation  levels  promote  an 
increase  in  the  power  output  and  of  the  drying  efficiency, 
respectively. 

-  The  exergetic  losses  were  very  significant,  mainly  because  the 
hot  airflow  generated  was  not  used  to  drive  a  turbine,  therefore 
there  is  a  great  potential  of  exergy  to  be  harnessed. 
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